The Adsorption and Collapse Transitions in a Linear Polymer Chain near an 

Attractive Wall 
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We deduce the qualitative phase diagram of a long flexible neutral polymer chain immersed in 
a poor solvent near an attracting surface using phenomenological arguments. The actual positions 
of the phase boundaries are estimated numerically from series expansion up to 19 sites of a self- 
attracting self avoiding walk in three dimensions. In two dimensions, we calculate analytically phase 
boundaries in some cases for a partially directed model. Both the numerical as well as analytical 
results corroborate the proposed qualitative phase diagram. 
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I. INTRODUCTION 



The behaviour of flexible polymers in solution at large 
length scales is independent of the chemical nature of 
the polymer and the solvent, and these universal scal- 
ing properties are well understood in terms of the renor- 
malisation group approach [0,0 . The polymer chain is 
known to undergo a transition from a random-coil phase 
to a globular phase as the temperature, or the pH of the 
solution is varied. The model of a self-avoiding walk on 
a lattice with nearest-neighbour attractive or repulsive 
interactions provides a simple model to understand the 
collapse transition in polymers 0. 

When the chain interacts with an impenetrable surface 
its conformational properties are strongly modified [|4|,|3| . 
There is a competition between the lowering of internal 
energy near an attractive wall, and the loss of entropy 
due to constraints imposed by the impenetrable surface. 
For a strongly attractive surface, the polymer sticks to 
the surface, and for weak attraction, the preferred state 
is away from the surface. Thus, there is a transition 
from the state when the chain is mostly attached to the 
surface to a detached state, when the temperature is in- 
creased. This behaviour finds applications in lubrication, 
adhesion, surface protection, etc M . 

If there is self-attraction as well as attraction to the 
surface, there is the possibility of a collapse transition 
in the desorbed polymer, or in the polymer adsorbed on 
the surface. In addition, there is a surface-attached glob- 
ular (SAG) phase, in which the polymeric globule gets 
attached to the attractive surface J7|. In the thermody- 
namic limit, the SAG phase has the same free energy per 
monomer as the bulk globular phase, and the transition 
between them is a surface transition. In earlier papers 
[Mm, we had discussed the phase diagram in this case, 
and investigated the phase diagram in a lattice model us- 
ing extrapolation of exact series expansions. This scheme 
has been found to give satisfactory results as it can take 
into account the corrections to scaling. To achieve the 



same accuracy by the Monte Carlo method, a chain of 
about two orders of magnitude longer than in the exact 
enumeration method has to be considered J9[. 

In this paper, we show that the qualitative features of 
the phase diagram in three dimensions can be determined 
by simple phenomenological arguments. In the case of a 
partially directed polymer in two dimensions, we deter- 
mine the exact phase diagram of the SAG phase. In this 
case, the polymer has different behaviour depending on 
whether it is near the wall perpendicular to the preferred 
direction (SAG1) or the wall parallel to the preferred di- 
rection (SAG2). We determine the phase boundaries of 
SAG1 and SAG2 phases by calculating their orientation 
dependent surface energy. We also determine the tran- 
sition between SAG1 and SAG2 phases when both walls 
are present. We also summarise our results of analysis 
of exact series expansion in three dimensions which we 
have extended by two more terms. 

The paper is organised as follows. Sec. contains the 
definition of the model and of the various phases. In 



Sec. Ill, we briefly review earlier work before providing 



arguments for the qualitative nature of the phase dia- 
gram in two and three dimensions. The phase diagram 
obtained is compared with numerical results from series 
expansion in Sec. IV. Sec. M contains the analytical re- 
sults obtained for the partially directed model. 



II. MODEL AND DEFINITIONS 



A simple lattice model for a linear polymer in a poor 
solvent is a self avoiding walk (SAW) on a regular lattice 
with an attractive interaction energy e„ between pairs of 
sites of the walk which are unit distance apart but not 
joined by an step of the walk. The adsorbing surface is 
modelled by restricting the walk to lie in a upper half 
plane and by associating an attractive energy e s with 
each monomer [site of the walk] lying on the surface. 
In the partially directed self-avoiding walk (PDSAW) in 



two dimensions, there is an additional restriction that the 
walk cannot take steps in the negative x-direction. 

We will work with the reduced variables ui = e@ es and 
u = e^ 6 ", where j3 is the inverse temperature. For the 
sake of easy reference, we now define all the phases that 
we will encounter later on, at one place. Consider a poly- 
mer consisting on N monomers whose one end is attached 
to a fixed site on the surface. If e s and e u are small in 
magnitude, the polymer exists in the swollen random-coil 
phase, away from the surface. In this phase, the mean 
radius of gyration varies as N v where v takes the self 
avoiding walk value \v w 0.588 in 3D and v = 3/4 in 
2D]. The number of monomers in contact with the sur- 
face is of order one in this case. We shall call this phase 
the desorbed extended (DE) phase. If e u is large and 
e s is small, the polymer exists away from the wall as a 
compact ball of finite density. In this case, the radius of 
gyration of the polymer varies as JV 1 ' in d-dimensions. 
We shall call this phase the desorbed collapsed (DC) 
phase. If the surface attraction e s is sufficiently large, 
the polymer sticks close to the surface. In this case, a fi- 
nite fraction of monomers are on the surface, and the ex- 
tent of the polymer perpendicular to the surface is finite. 
Along the surface the polymer roughly acts as a polymer 
chain in (d — l)-dimensions. Depending on whether the 
attractive self-interaction is large or small, the polymer 
is in a collapsed phase, with its transverse size varying 
as N 1 ^ d ^ 1 \ or in the extended phase with the trans- 
verse size varying as N v , where v' is the self avoiding 
walk exponent in (d — l)-dimensions. We shall call these 
phases the adsorbed collapsed (AC) and the adsorbed 
extended (AE) phases respectively. In addition to these 
phases, the polymer may exist as a collapsed globule of 
finite density which sticks to the surface. In this case, 
the size of the polymer in the directions transverse and 
perpendicular to the surface varies as N x / d and the num- 
ber of monomers in contact with the surface varies as 
jy(d-i)/d_ ^y e g^gjj c& \\ t ms phase the surface adsorbed 

globule (SAG) phase. Note that in 2-dimensions, the 
AC and the AE phases cannot be distinguished from one 
another. 

The polymer undergoes a transition from the extended 
to the collapsed phase when the temperature is varied. 
At the transition temperature between the DC and the 
DE phases, called the 0-point, the critical behaviour is 
described by a tricritical point of the O(n) (n — » 0) spin 
system. At the 0-point, R b ~ N ve with vg = 4/7 for 2D 
(U| and 1/2 for 3D 0. The transition from AE to AC is 
described by vg corresponding to one lower dimension. 



III. QUALITATIVE PHASE DIAGRAM 

We first briefly review earlier results obtained for the 
model when both the monomer- monomer interaction as 
well as the interaction with the wall are attractive. One 



of the earlier papers on such systems was by Bouchaud 
and Vannimenus |ljj , in which the the phase diagram 
of a polymer living on a Sierpinski gasket was analyti- 
cally derived. The phase diagram consisted of the AE, 
DE and the DC phases. In |12[ , the phase diagram in 
2-dimensions was obtained approximately by series ex- 
pansions and it was found to be qualitatively similar to 
that for the gasket. In M, the possibility of the existence 
of the SAG phase in 2-dimensions was pointed out based 
on analysis of series expansions. Evidence for the exis- 
tence of a surface transition from the SAG to DC phase 
was also presented. A variant of the model, the PDSAW 
model in 2-dimensions, has been more amenable to an- 
alytical calculations. For PDSAW, in 2-dimensions, the 
exact calculation of the phase boundary between the col- 
lapsed and the extended phases []13| — |l5|| has been numeri- 
cally confirmed in |L6| . The phase diagram thus obtained 
is qualitatively similar to the undirected 2-dimensional 
model. In |L7J, the existence of the SAG phase in the 
PDSAW has been argued for, based on series expansion 
analysis. 

The model is less studied in 3-dimensions. Monte Carlo 
simulations |18| and series expansion analysis M on the 
cubic lattice showed the existence of four phases: AE, 
AC, DE and DC. While Q claimed the existence of two 
multicritical points, the earlier preliminary results B ob- 
tained from series expansion seemed to support one mul- 
ticritical point. More careful analysis of the series, re- 
ported later in the paper, shows that there are indeed 
two multicritical points. The question of whether an 
SAG phase exists in 3-dimensions or not has not been 
addressed so far. Also, the possibility of surface transi- 
tions among the collapsed phases has not been explicitly 
dealt with. Thus, in spite of many earlier studies, the 
qualitative behaviour of the system is not fully estab- 
lished. 

We now determine the qualitative nature of the phase 
diagram from phenomcnological considerations. Con- 
sider the case when the interaction with the wall is repul- 
sive, i.e., u) < 1. The polymer will be in a desorbed state 
because proximity to the surface results in increase of in- 
ternal energy as well as loss of entropy. As u is increased 
from 1 to oo, the polymer undergoes a collapse transi- 
tion from a desorbed expanded state (DE) to a desorbed 
collapsed state (DC) at a critical value u^ d (see Fig. [l]). 
When the interaction with the wall is attractive, there is 
a competition between lowering of free energy and loss of 
entropy as the wall is approached. For e s > 0, the loss of 
entropy is more dominant and hence the polymer remains 
desorbed. Now, as u is varied, the polymer undergoes a 
transition from DE to DC at the same critical value u^ d . 
Therefore, in the lower part of the phase diagram, there 
is a vertical phase boundary u — u c separating the DE 
and the DC phases. 

Consider now the case when the interaction between 
nearest neighbour monomers is close to zero, i.e., u ^ 1. 



Clearly, the polymer will be in the DE phase. Now, as lu 
is varied from 1 to oo the polymer undergoes a transition 
from DE to one in which it is adsorbed onto the wall and 
extended, i.e., the AE phase. Let this transition occur at 
a critical curve lu c (u) that intersects the w-axis at lu*. 




FIG. 1. The qualitative phase diagram in three dimensions. 



Now, consider possible phases of the system for large w 
and u. At T = 0, the polymer has the maximum possible 
density. It can then be modelled as a Hamiltonian walk. 
The attractive energy per site is — (<i — l)e u . In addition, 
for a finite polymer of N monomers, there is a surface 
energy cost which is easily seen to be d€ u N t ^ d ^ 1 ^ d . Then 
the free energy of the DC phase at T = is 



E DC = -( d - l ) e " N + d ^N^ 



(1) 



At T = 0, for lu c i < lu < cu c2 (see Fig. |j) the poly- 
mer exists in the SAG phase as a rectangular paral- 
lelepiped of size L\\ and L± in directions parallel and 
perpendicular to the surface. Its bulk energy is the 
same as in the DC phase and the surface energy is 
(e u — e^Ljj^ 1 + (d — l)e u N/L\\. Minimising the surface 
energy with respect to Lu, we obtain 



V SAG = -(d-l)e u N + de u d {tu-e s )-*N- 



(2) 



In the AC phase, we have L±_ = 1, L\\ = iV 1 /^-!) an d 
the free energy at T — is 



E 



AC 



-{d-2)e u N-e s N+{d-l)e u N—. (3) 



From Eqs. fl] and @, it follows that the surface energies 
become equal when e s = 0, i.e., lu c \ — > 1 as u — > oo. 
Also, the bulk free energy term for the AC phase has 
lower value than that of the SAG phase when e s > e u . 
Thus, lu c2 — > u, as u — > oo. 

Consider now the case when lu = oo. Then the poly- 
mer is adsorbed onto the (d— l)-dimensional surface. In 
d > 2, there is the additional possibility of the AC phase. 



As u is increased from 1 to oo, the polymer will undergo a 
collapse transition, while still remaining adsorbed, from 
AE to AC. This transition occurs at the critical value 
of (d — l)-dimensional collapse u* d _ l . Since the critical 
temperature for collapse transition increases with dimen- 
sion, u* l _ l > u* d . We will now argue that the curve lu c3 
originating from (m^, °°) bends to the right when lu is 
decreased to a finite but large value. The partition func- 
tion, when written as perturbation series in lu^ 1 , is 



Z(u, lu) — Zq(u)lu 



N 



N 
~1 



n 



Hi 
u 



«2 
,,2 



(4) 



where rij is the fraction of bonds such that their end 
points have exactly j nearest neighbour monomers. 
Clearly, no is larger in the AE phase as compared to 
the AC phase, while ri\ and n^ are smaller. Using 
tiq = 1 — m — n2 in Eq. EL it follows that for large but 
finite lu, the free energy is lower for the AE phase. Hence, 
the phase boundary lu c3 curves to the right. 

The phase diagram for the 2-dimensional problem is 
qualitatively the same as that of the 3-dimensional prob- 
lem except that there is no AC phase, and hence no lu c3 
phase boundary. We now argue that the phase bound- 
aries u Cl lu c , lu c \ and lu C 2 meet at one point. For the 
sake of clarity, we will illustrate the arguments for the 
2-dimensional problem. In the DC and the SAG phases, 
the polymer is a compact two dimensional object with 
finite density. We define a (8) as the surface tension be- 
tween the surface of this object and the liquid, where 8 
is the angle the surface makes with the horizontal. For 
a shape r(8), the free energy is a sum of two terms: the 
bulk term which depends on u alone and a surface term, 
which can be written as an integral over the angle depen- 
dent a (8). 

Near the phase boundary lu C 2 separating the AE and 
the SAG phases, the shape is highly anisotropic and 
R s ^$> Rb, where R s and Rb are the extent of the polymer 
along and perpendicular to the surface. R s diverges as 
we approach the phase boundary from within the SAG 
phase. Additional cost of creating two surfaces of ori- 
entation 8 — should be zero. Hence, along the phase 
boundary lu C 2, we have 



cr(0) + a w = 0, 



(5) 



where a w is the free energy cost per unit length when 
the polymer is along the wall. Near the phase boundary 
lu c \ separating the DC and the SAG phases, the shape of 
SAG is such that the part in contact with the wall has 
orientation 8 = 0. Clearly, this configuration becomes 
unfavourable in comparison to the DC phase when 



cr(0) = cr u 



(6) 



For the DE — DC transition, clearly the surface tension 
must vanish at the collapse point. Thus, along u c we 
have 



o-(O) =0. 



(7) 



From Eqs. |5|-[7|, it is clear that if any two of the phase 
boundaries intersect at some point, then the third will 
also pass through that point. 

It still remains to be shown that uj c will also pass 
through the same point as the other phase boundaries. 
Let u and u> get transformed to u' and u/ under a scale 
transformation as 



u' = /(«), 



■jj 



(8) 
(9) 



The function f(u) is independent of the surface param- 
eter uj because u is a bulk parameter. There will be 
three fixed points for Eq. |8| namely u = 0, u = u* and 
u = oo where u* is the only repulsive fixed point. Con- 
sider Eq. y when u is fixed at each of its three fixed points. 
In the simplest scenario, for each value of u, there will 
be three fixed points of Eq. |9j, one corresponding to no 
attraction, one to very strong attraction and the third 
a repulsive fixed point. After fixing the flow directions, 
the final flow diagram looks schematically as shown in 
Fig. g. The fixed points Ai : A2, C\ and C2 correspond to 
the four phases. The fixed points A,Bi,C and Bi cor- 
respond to the four critical phases corresponding to the 
phase boundaries and the point B corresponds to the 
repulsive multicritical point. 



at the origin and all other sites are confined to the half 
plane y > 0. Let Cn(N s ,N u ) be the number of SAWs 
of N sites having N s monomers on y = and N u near- 
est neighbour monomer pairs. In jq], we reported the 
enumeration and analysis of the series Cn(N s , N u ) up to 
N = 17 for the cubic lattice. We have now extended the 
series for 3-dimensions by two terms and reanalysed the 
data to obtain a better estimate of the phase boundaries. 

For fixed u, we identify the position of the phase 
boundary separating the desorbed phase from the ad- 
sorbed or attached phases as that value of u> at which 

i *' is a maximum. Fig. shows the variation of 1 
for two values of u for N = 19. 
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FIG. 3. The dependence of 1 "' on uj is shown. For 
u = 2.0, there is only one peak corresponding to the DE to AE 
transition. For u = 3.5, there are two peaks corresponding to 
the DC to SAG to AC transition. 




FIG. 2. The schematic flow diagram. 



For fixed uj, we identify the position of the phase 
boundary separating the extended phase from the col- 



lapsed phase as that value of u at which 

imum. Fig. ^ shows the variation of 
of uj for N = 19. 



9{N. 



is a max- 



9{Nu) 



for two values 




IV. SERIES EXPANSION RESULTS 

We enumerated all SAWs up to a certain length on 
the cubic lattice in which the first site of the walk lies 



FIG. 4. The dependence of 1 "' on u is shown. 



For 



u — 2.0 there is only one peak corresponding to the DE to DC 
transition. For uj — 3.8, there are two peaks corresponding to 
the AE to AC to SAG transition. 



The values of u^ d and uj* obtained by this method is 
2.00 and 1.49 respectively. The previous results were 
1.76 by series expansion method H and uj* — 1.45 



*3d 



by Monte Carlo method |19| and uj* — 1.5 by series ex- 
pansion method [E0|. It is possible to obtain better esti- 
mates of Mg d as well as the phase boundaries by extrap- 
olating for large N. Let 



Z n (uj,u) 



E 

N S ,N V 



C n (N s ,N u )lj 



N, U N„ 



At zero temperature, it is easy to see that the configura- 
tional energy of the polymer is minimised if it assumes a 
square shape of size y/N by \N. For small nonzero tem- 
peratures, the polymer assumes a shape which is slightly 
perturbed from this zero temperature square shape. We 
will derive an effective surface energy for these fluctua- 
tions in Sec. VA. Using these results, we determine the 
shape of SAG1 and SAG2 in Sec. |VB|. In Sec. fV(l we 



(10) calculate the phase boundary between the various phases. 



be the partition function. Then, the reduced free energy 
per monomer can be written as 



1 



G{u) 1 u)= lim — log Z n (uj,u) 



N 



N 



(11) 



We refer to MM for details of the methods used for ex- 
trapolating to large N in Eq. nTl The phase boundaries 



arc then found from the maxima of 



G(ui,u) ,_ d(N s ) 



<1<- 



d<L s 



and 5^*1 ( 



d(N u ) ■ 



The phase diagram thus obtained is shown in Fig. H. 
We obtain u^ d = 1.76 and uj* — 1.48 which accords fairly 
well with the previously obtained results. The phase di- 
agram obtained from series analysis agrees qualitatively 



with the phase diagram proposed in Sec. Ill 
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FIG. 5. The phase diagram for a polymer in 3-dimensions 
obtained from series analysis is shown. The dashed and the full 
lines are the locations of maxima of — „ ^ and — _ ■ ■ , re- 
spectively. 



V. ANALYTIC CALCULATION FOR THE TWO 
DIMENSIONAL DIRECTED POLYMER 

In this section, we determine the phase boundary sepa- 
rating the SAG phase from the DC and the AE phases in 
the PDSAW model. We do so by calculating the macro- 
scopic shape of the collapsed phases at low temperatures. 



A. Effective surface energy 

For the directed polymer in the collapsed or the SAG 
phases, the density in the bulk is exactly one and the con- 
figuration is "frozen" . Only the position of the boundary 
can change, as there is some fluctuation of height allowed 
at the boundary. Thus /^^r, (SAG) = /bulk^^) = — e «> 
independent of uj. Consider a polymer shape as shown in 
Fig. |. 
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FIG. 6. Schematic diagram of a partially directed polymer 
for T > 



The energy of the configuration is 



b-2 



E = -e u N + f(ai + a 2 + 26) + -f "£ \y J+2 - 



Vj 



3=0 



(12) 



By a redefinition of E, we drop the bulk term propor- 
tional to N. The shape of the polymer is determined 
by the rest of the terms which are all proportional to 
vJV. We replace the terms under the summation by an 
integral over an effective orientation dependent surface 
energy /(#), where 9 is the angle the surface makes with 
the horizontal. In this case, it is straightforward to cal- 
culate f(9). Consider all possible walks with an average 
slope tan(#) = y/x. Then, the sum over all weighted 
paths is 



-/3xsec(9)f(6) _ 



E ME*-* IL 



(13) 



Vl,--;Vx 



where p = e - ^ and 5 is the usual Kronecker delta func- 
tion. Taking Laplace transform with respect to y, we ob- 
tain independent summations over j/^. These are easily 
done giving 



m - \ 



sin6»log(z ) 



cos 9 , (z 

—FT- l0 § 



p)(l-pz ) 



zo(l-p 2 ) 



where 



zq = 



(1 



- p 2 ) tan 6 + v/(l-P 2 ) 2 tan 2 
p(l + 2tanfl) 



. p^ 



(14) 



(15) 



We also need to calculate the energy cost a w of ad- 
sorbing onto the wall unit length of the polymer. For 
SAG1, it is trivially equal to a w \ — —e s . We calculate 
a W 2 for SAG2 by the transfer matrix method. If ipi de- 
notes the y-coordinate of the lowest portion of the poly- 
mer at site i, then the weight of obtaining ipi+2 from ipi 
is (^i|T|Vi +2 ) = [l + (^ 2 -l) ( 5. 0i+2 , o ] U -l^-^+-l/ 2 . By 
trying out an ansatz tpi = a + Sifiipo for the eigenfunc- 
tion, it is not difficult to verify that the largest eigenvalue 
of the transfer matrix T is 



A 



u 2 {u 2 -l)(u-l) 

LU 2 (u — 1) — U 



Then, a w2 — — log(A)/(2/3). Clearly, as u 
the correct limit — e s . 



(16) 
co, a W 2 has 



C. Phase Diagram 



We calculate the phase diagram for the directed poly- 
mer from Eqs. g-Q. These equations give most of the 
phase boundaries except the transition between SAG1 
and DC and the transition between SAG1 and SAG2. 
The SAG1-DC transition is not given by Eq. || because 
the shape in contact with the surface does not have orien- 
tation 8 — 0. This anomaly arises due to the constraint 
of directedness. The surface transition from SAG1 to 
SAG2 is one in which the globule would have lower free 
energy if attached to the x- wall rather than the y-wall. 




FIG. 7. The shape of SAG polymer is shown for different 
values of u) when u is kept fixed at 10.0. The position of the 
wall is denoted by a dotted line (vertical for SAG1 and hor- 
izontal for SAG2). The shape of SAG1 corresponds to the 
part of the curve from the wall to the right, while the shape 
of SAG2 corresponds to part of the curve above the wall. 



B. Calculation of the macroscopic shape 



In this subsection we describe the shape determined 
by minimising the surface energy of the collapsed phases. 
Given the expression for the temperature and orientation 
dependent f(8), and also the value of surface energy of 
polymer attached to the wall, it is straightforward to de- 
termine the globular shape which minimises the surface 
energy given a fixed volume. This is the classical Wulff 
construction. The result is that the macroscopic shape 
of the polymer is given by 



2/3Ay 



c 2 e pAa; (l-pc 1 e- pAa; )(c 1 e 



/3\x\ 



-/3\x 



■P), (17) 



where the two constants c\ and ci are fixed by the two 
boundary conditions. The Lagrange multiplier A is de- 
termined by the constraint that the total area under the 
curves is N. The constants C\ and c 2 are now varied to 
obtain the shape with the lowest surface energy. 

Fig. M shows the shape of the SAG polymer for differ- 
ent values of u>. All the shapes lie on top of each other if 
we scale the coordinates as X = (3\x and Y — j3\y. 



Transition from DC to DE (u c ): The critical value 
u c is obtained from Eq. [7J i.e, cr(0) = 0. This is equiv- 
alent to the 2/(0) + e u = 0. Substituting for /(0), we 
obtain 



y^-1 



1 

j 
u 



which has the solution 



3.38298. 



(18) 



(19) 



Note that this result matches exactly with the result 
for DC-DE transition obtained by the transfer matrix 
method |TJ-jl§. 

Transition from SAG1 to AE (w C 2) s This phase 
boundary is determined by equating the perpendicular 
extent of the SAG1 phase above the wall, to zero. This 
condition gives rise to the phase boundary 



Wc2 



1 + u 2 + y/(l + u 2 ) 2 - 4u 3 



2u 



(20) 



This solution has a natural boundary at u = u c at which 
value the expression under the square root sign becomes 
equal to zero. 

Transition from SAG1 to DC (oj c i) '■ The transi- 
tion from SAG1 to DC occurs when the energy cost of 
creating a globule sticking to the wall becomes equal to 
the energy of a DC polymer. We omit the algebra and 
state the final result 



LU cl (u) 



1 



va 



4u 3 



2a 



(21) 



Previous analytical studies on the PDSAW |L3|-|15| had 
considered the case when the wall was only along the x- 
direction. The results obtained above for SAG1 are for 
a wall along the y-direction. While the numerical values 
for the phase boundary differ, the phase diagrams are 
qualitatively similar. 

Transition from SAG2 to AE (u c2 ) : From Eq. |, 
the phase boundary uj C 2 is given by a(0) + a w = 0. Sub- 
stituting the values of the surface energies and solving 
for bj. we obtain 



diagram obtained is qualitatively similar to the phase di- 
agram proposed in Sec. III. The additional transition 
between the SAG's is a consequence of the directed na- 
ture of the PDSAW model. 
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FIG. 8. The phase diagram for the 2-dimensional PDSAW 
model. The DE-AE phase boundary is schematic. 
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+ V7 



4u 3 
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(22) 



where a = 1 + ^/u — u 2 + u 5 ' 2 . The phase boundary 
u) C 2 has a natural boundary at u — u Cl at which value 
the expression under the square root sign becomes equal 
to zero. The result differs f rom the tra nsfer matrix re- 
sult H-f5|, w<a = *±i + ^ /(ll 2(t 1 -i)' 4 " 3 ■ However, this 
discrepancy is solely due to the fact that we consider 
only wall, while the transfer matrix approach required 
two parallel walls. This corresponds to changing Eq. |5| 
to 2a w2 + e u = 0. 

Transition from SAG2 to DC (w c i) : From Eq. |, 
this transition occurs when cr(0) = a w . The resulting 
equation can be solved to obtain 



(23) 



1 



This covers all the transitions when we consider SAG1 
and SAG2 separately. But if we consider the scenario 
where the possibility of both SAG's are allowed, then 
there is a surface transition from one to the other when 
u and w are varied. 

Transition from SAG1 to SAG2 : This transition 
is determined by equating the surface energies of SAG1 
and SAG2. However, it turns out that we cannot obtain 
a closed form expression for the phase boundary. Instead, 
we solved for it numerically using MATHEMATICA. 

In Fig. @, we plot the phase diagram when both SAG1 
and SAG2 are allowed to exist. Note that the phase 
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